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Abstract 

Renormalization group (RG) smoothing is employed on the lattice to investi- 
gate and to compare the monopole structure of the SU (2) vacuum as seen in 
different gauges (maximally Abelian (MAG), Polyakov loop (PG) and Lapla- 
cian gauge (LG)). Physically relevant types of monopoles (LG and MAG) are 
distinguished by their behavior near the deconfining phase transition. For 
the LG, Abelian projection reproduces well the gauge independent monopole 
structure encoded in an auxiliary Higgs field. Density and localization proper- 
ties of monopoles, their non- Abelian action and topological charge are studied. 
Results are presented confirming the Abelian dominance with respect to the 
non-perturbative static potential for all gauges considered. 
PACS: ll.15.Ha, 12.38.Gc 
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I. INTRODUCTION 



Over the last two decades various attempts have been aiming at a quahtative understand- 
ing and modehng of two basic properties of QCD: quark confinement and chiral symmetry 
breaking. Two well-known, complementary schemes are the instanton liquid model and 
the dual superconductor picture of the QCD vacuum While the first model explains 

chiral symmetry breaking and solves the Ua{^) problem, the second provides the language 
to discuss and to quantitatively describe the confinement mechanism within an effective field 
theory. Nowadays, it has become of central importance in the lattice community to quantify 
the properties of the lattice vacuum guided by the instanton liquid picture 0]. The evidence 
for a particular instanton structure is still controversial, in pure Yang-Mills theory and in 
full QCD. 

In the effective dual Abelian Higgs theory ||^, the vacuum is viewed as a dual supercon- 
ductor, where condensation of color magnetic monopoles (described by the Higgs phase 
of that model) leads to confinement of color charges by flux tubes through a dual Meissner 
effect. In order to establish this correspondence step by step, it has first been demonstrated 
that the Abelian components of the gauge fields are effective degrees of freedom (among the 
original non- Abelian gauge fields) at large distances, leading to the concept of Abelian dom- 
inance The role of monopoles and their different percolation properties for confinement 
and deconfinement were empirically substantiated by a large number of lattice simulations 
over the last years. It was shown that in the confinement phase monopoles percolate through 
the 4D volume |^ and are responsible for the dominant contribution to the string tension 
(monopole dominance) More recently, these monopoles have been shown to be pre- 

dominantly selfdual (in confinement) [jll| and to encode the non-perturbative structure of 



the gauge field at large distances as seen by quarks, explaining the chiral condensate |12 



as well as the topological density [T^. In order to establish the link to the effective dual 
theory, the action that describes the first quantized theory of monopoles has been derived 
from the configurations of magnetic currents |T^. At the same time, inherited from the 
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original gauge fields, these "Abelian monopoles" are characterized by a self-action (which is 
essentially no n- Abelian!) and a complicated interaction with the remaining degrees of free- 
dom. In the last couple of years also a deeper connection between monopoles and instantons 
has been pointed out, both on the lattice p5|-[T7[] and in the continuum [p!8| , pO| -p^,pJ[1 . 



Following 't Hooft [0], monopoles should be searched for as pointlike singularities of 
some gauge transformation dictated by a local, gauge covariant composite field. The choice 
of this field implies the choice of a particular gauge to analyze the non-Abelian gauge field. 
The best known examples are the maximally Abelian gauge (MAG) and the Polyakov 
gauge (PG). It has become practice to identify the monopoles as particle-like singularities 
on the lattice after Abelian projection from one of the gauges chosen, identifying their 
trajectories in the resulting compact Abelian lattice field. This is known as the DeGrand- 



Toussaint (DGT) construction [|2^. An attractive alternative, related to the Laplacian 
gauge (LG), has been revived recently [^,^. This method offers the possibihty to identify 
monopoles on the lattice without the problems that the conventional methods have and 
bears a close relationship to the 't Hooft-Polyakov monopoles. One can define them without 
really performing a gauge transformation. The LG can, however, also be the starting point 
for an Abelian projection followed by a DGT monopole construction. In this paper we intend 
to compare these three gauges used to define DGT monopoles with respect to their physical 
significance at the deconfinement transition. Furthermore we want to give a justification of 
the DGT construction of monopoles in the case of the LG, in view of their gauge independent 
localization similar to that of 't Hooft-Polyakov monopoles. 

An important ingredient of our discussion is a smoothing method for non-Abelian gauge 
fields used to remove UV fluctuations that are without relevance for the long distance physics. 
From our previous experience with this method we know that this procedure is suitable to 
eliminate UV lattice artifacts from the monopole-related observables, too, without destroying 
the confining structure as usual cooling procedures finally do. Originally, methods like 
cooling or smoothing have been developed in order to study the instanton structure which 
otherwise is possible only by fermionic (spectral flow) methods |2^. A particular feature 



of cooling or smoothing, that should interest us here, is that it maps gauge fields on full 
gauge field configurations which still can then be put into various gauges for a closer study 
of one or the other qualitative picture. Concerning the monopole degrees of freedom, it 
is advantageous that smoothed configurations carry monopole currents with short loops 
removed. 

The renormalization group (RG) based smoothing method which uses (classically) perfect 



actions has been strongly recommended in the last two years because it is not 



afflicted with some problems of the simple or improved cooling methods. UV fluctuations 
near to the lattice discretization scale are removed from the gauge fleld by one blocking step 
followed by a smooth interpolation of the fleld on the original lattice. By this procedure, 
performed locally below a well-deflned (the blocked) scale, the RG smoothed conflgurations 
remain conflning (if they are derived from MC runs in the conflning phase). The string 
tension is reproduced almost completely. At the lattice sizes (and values of P) that were 
to our disposal, the two RG levels involved in the procedure are not totally decoupled from 
the conflnement scale. Fluctuations are removed by the blocking-inverse-blocking scheme, 
which contribute a few percent to the string tensions at the distances where we can study 
it. 

In this paper we put special emphasis on the possibility of a gauge invariant deflnition 
of monopole trajectories and test it against the conventional DGT method usually applied 
to conflgurations being Abelian projected from various gauges. We were particularly inter- 
ested in temperatures below and above the deconflnement critical point. This gave us the 
opportunity to compare the behavior of the respective monopole degrees of freedom at the 
transition and examine their dynamical relevance. 

The paper is organized as follows. In Sec. II we briefly discuss the RG smoothing tech- 
nique and recall some properties of action and topological charge densities after smoothing. 
In Sec. Ill we describe the different gauge flxing procedures being used and report certain 
technical details how they work for smoothed conflgurations. In Sec. IV various aspects 
of the monopole degrees of freedom are presented, pinning down the main difference be- 



tween the different gauges by means of the different physical behavior of the corresponding 
monopoles across the phase transition. In Sec. V we present some incidental results of our 
study concerning Abelian dominance, comparing the Abelian string tensions obtained by 
Abelian projection from various gauges with the fully non- Abelian non-perturbative poten- 
tial extracted by smoothing. In Sec. VI we summarize our conclusions. 



II. RG SMOOTHING AND PHYSICAL DENSITIES 

To recognize semiclassical structure in ensembles of gauge field configurations generated 
by lattice simulations, cooling methods have been used, that minimize the action. However, 



even improved versions of cooling |]31[ gradually reduce the string tension. These versions 
are designed to stabilize instantons above a certain threshold for their scale size p (typically 
p > 2a with the lattice spacing a). However, they let close instanton-antiinstanton pairs 
annihilate. 

Up to now, most lattice studies are performed using the Wilson action, even if they 
focus on the topological structure of the vacuum. This action is known to be afflicted with 
so-called dislocations, structures of a size near one lattice spacing. For dislocations certain 
lattice definitions of the topological charge Q give a signal of unit charge, although their total 
action 5* violates the inequality S > \Q\87i'^/g'^ and the entropic bound S > S = ASn^/ll N^. 
For discretized instantons, the Wilson action decreases with size p of an instanton, such that 
isolated instantons shrink under cooling before they decay as a dislocation. Thus, cooling 
with Wilson action is unsuitable to unambiguously determine the instanton size. Other 
methods like APE smearing let instantons even grow. 

To avoid these ambiguities we have proposed to use the renormalization group motivated 



method of "constrained smoothing" which is based on the concept of perfect actions p2 
For discretized instantons, these actions guarantee a size independent action for instantons, 
even with radii near to the lattice spacing. They allow a theoretically consistent "inverse 
blocking" operation which, for instance, reconstructs an instanton solution on a finer lattice 
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with the same value of action. Inverse blocking is a method to find a smooth interpolating 
field on a fine lattice by constrained minimization of the perfect action, provided the config- 
uration is given on a coarse lattice. This way to interpolate allows to define an unambiguous 
topological charge p^ . 

For non-classical configurations like generic Monte Carlo configurations, one is interested 
to study smoothed configurations which are made free from UV fiuctuations below a well- 
defined scale. Constrained smoothing in the way we use it, consists in first blocking the fields 
{U}, sampled on a fine lattice with lattice spacing a using the perfect action, to a coarse 
lattice configuration {V} with lattice spacing 2 a by a standard blockspin transformation. 
Then inverse blocking is used to find a smoothed field {W™} replacing {U} on the fine 
lattice. We refer to this method as "renormalization group" (RG) smoothing. 

This procedure can be cyclically repeated, choosing different blocked lattices for the 
blocking step in subsequent cycles. This has been practized in Ref. |3^. Experience shows 
that this RG cycling method makes configurations locally more and more classical. Unfor- 
tunately, there is no well-defined smoothing scale, but in contrast to cooling this method 
preserves long range features like confinement rather well. An important advantage of the 
simpler RG smoothing method outlined above is that it does not drive configurations into 
classical fields as unconstrained minimization of the action (and RG cycling) would do. It 
saves the long-range structure of the Monte Carlo configuration in {V} and allows to study 
semiclassical objects in the smoothed background deformed by classical and quantum inter- 
action. It is the upper blocking scale which roughly defines the border line between "long" 
and "short range". 



In this work we used a simplified fixed-point action p9|j30[| for Monte Carlo sampling and 
for RG smoothing. Before describing in detail the various gauge conditions we should explain 
which gauge invariant quantities we shall use to characterize the smoothed configurations. 
Partly they are suggested by the perfect action Spp itself which is parametrized in terms of 
only two types of Wilson loops, plaquettes Uci = Ux,^jt,v (type Ci) and tilted 3-dimensional 
6-link loops (type C2) of the form 
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and contains several powers of the linear action terms corresponding to each loop of both 
types that can be drawn on the lattice 

Sfp{U) = E EE^(^,J')(1 - ^ Tr . (2) 

type i d j=l 

The parameters of this action are reproduced in Table |. 

We define the action density Ssite{x) per lattice point (a local action) as follows 

ssueix) = E f E - I Tr Uc,^^,y + E - I Tr Uc^.^A . (3) 

i=l \Ci(x) C2(x) / 

Here, Cj{x) {j = 1,2) means loops of type j running through the lattice site x. Summing 
Ssite{x) with respect to x yields the total action of the configuration. 

The topological charge density definition we are using is (i) the simplest "field theoretic" 
one constructed out of plaquettes around a site x 

1 +^ 

^(^) = E ^f^uaptr {U^,f,,Mx,a,p) (4) 

^ p,,i^,a,p=-4: 

and (m) the charge contributed by hypercubes according to Liischer's definition of topological 



charge p3. We did not attempt to improve the field theoretic definition. Both topological 



densities used are known to behave regularly for smooth configurations |]S5 

In order to give a characterization of Abelian monopole currents in terms of locally 
defined gauge independent observables, we consider also the 3-cubes of the original lattice 
dual to an elementary piece of monopole world line (carried by a link of the dual lattice). 
This leads us to a natural definition of a local action on the monopole world line, ss-cube{c) 
instead of SsUei^). We include into s^^cubeic) the contribution of all plaquettes forming the 6 
faces of the 3-cube c plus the contribution of all 6-link loops which wind around its surface. 
In short, it is the part of the total action which lives on that cube. 



III. MONOPOLE IDENTIFICATION METHODS 
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A. Maximum Abelian gauge and Laplacian gauge 



The most popular gauge is the maximally Abelian gauge (MAG) [^]. Abelian projection 
of configurations put into MAG exhibits Abelian dominance. This gauge is often used to 
study the dynamics of monopoles on the lattice. It is enforced by an iterative minimization 
procedure, which can get stuck in local minima that are gauge copies of each other (so-called 
technical Gribov copies). The Laplacian gauge (LG) is not afflicted with this problem. The 
basic idea has been suggested more than a decade ago but has apparently not been used 
until very recently [^]. Originally, in the context of the adjoint Higgs (Georgi-Glashow) 
model, it has been proposed to introduce, besides the quantized scalar matter field, a purely 
auxiliary adjoint slave field. 

Such a Higgs field can be used to define a gauge transformation to the "unitary" gauge 
with respect to the auxiliary field except at its zeroes. Thus, the corresponding monopoles 
are located at the zeroes of the auxiliary Higgs field and are representing the singularities 
of the gauge rotation. They can be identified as t'Hooft-Polyakov monopoles. For each 
given (generic, non-classical) non-Abelian gauge field the auxiliary field is determined as 
the eigenvector related to the lowest eigenvalue of the adjoint covariant Laplacian. This is 
why the resulting gauge is named Laplacian gauge (LG). The monopole identification via 
the Higgs zeroes is obviously gauge invariant (see also Ref. I^^). It is this prescription we 
want to apply to pure Yang-Mills theory on the lattice and to compare with other ways for 
the Abelian projection. Strictly speaking, in order to describe the location of the monopole 
world lines (Higgs zeros), we would need algorithms interpolating the Higgs field between 
the lattice sites. For our present purposes it will be sufficient to detect clusters of lattice 
points with small Higgs modulus. 

Technically, MAG and LG are closely related to each other. For MAG, the gauge func- 
tional to be minimized can be written as follows 
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with the gauge transformation VL^ acting on {f/} 



E - E ^X^i^ - bD.Xf , (5) 



r/(^) = Q U ~ 



encoded in an auxiliary adjoint Higgs field 

a 

subject to local constraints Y^ai^xY — ^ ^^"^ with adjoint links 



To change from MAG to LG, the local constraints ^ai^xY — 1 relaxed and replaced 
by a global normalization: Y^x^ai^xY = ^5 such that Eq. (|) can be further written: 

/ {D,xf - EE^:{-°";t(^)}^J • (6) 

Thus, the minimization of the gauge functional is reduced to a search for the lowest eigen- 
mode of the covariant lattice Laplacian. The LG of a given lattice configuration is un- 
ambiguously defined, except for the case that degenerate lowest eigenmodes exist. The 
corresponding LG configurations would then define true Gribov copies of each other. For 
both MAG and LG, the gauge transformation is finally accomplished by finding Vt^ that 
diagonalizes the field $3,. 

Here a technical observation ought to be made. As usual, the MAG is found iteratively 
by minimizing the functional (^. The fact that we analyze smoothed configurations does 
not mean that the MAG is computationally less demanding to find. Still, as in the case for 
unsmoothed configurations, we need 0(1000) iterations to fulfill our stopping criterion. We 
require that the next gauge transformation should deviate from unity uniformly by less than 
10^''. By abelianicity R we denote the average Abelian fraction per link after completion 
of the minimization. We show the (3 dependence of this quantity in Fig. |I] in comparison 
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with other gauges which do not exphcitly attempt to get a large abehanicity. There is 
a monotonous increase with (3 across the phase transition, but no "deconfinement signal" 
related to it. 

For the search of the lowest eigenvalue for the LG and its corresponding eigenvector, we 
apply the conjugate gradient algorithm to minimize the following Ritz functional 

_ (X, {-D^)X) 
^- {X,X) ' 

where is the adjoint lattice Laplacian. We choose 5 random start vectors = X^aa and 
follow the relaxation of the estimated eigenvalue over 3000 iterations. The vector with the 
lowest eigenvalue is then taken further representing the given gauge field configuration. The 
average over an ensemble of 100 gauge fields (per /? value) of the respective lowest eigenvalue 
is shown in Fig. |^ as a function of f3 (or temperature). This average is f3 dependent in both 
phases with a steep drop across the transition. The corresponding abehanicity in Fig. |l| is 
always smaller than the abehanicity for MAG, but approaches it with increasing (3 (in the 
deconfined phase). Among the different gauges, the biggest rise of abehanicity accompanying 
the transition is found in the LG case. 



B. Polyakov gauge 

This gauge condition puts special emphasis on the (gauge group valued) Polyakov loop 
P{x), the gauge transporter along the shortest loop starting from and arriving at a given 
lattice point which is closed by periodic thermal boundary conditions. The Polyakov gauge 
is enforced by simultaneous diagonalization of P{x) at each point. In Fig. ^ we also show 
the abehanicity achieved in the PG. Except near the transition (where it has a shallow 
maximum) the abehanicity is practically (3 (temperature) independent. This gauge has the 
smallest abehanicity among the discussed gauges, but it might be surprising that it comes 
near to 95 %. 
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C. Monopole identification 



After the gauge of choice has been fixed one extracts the Abehan degrees of freedom by 
Abehan projection. The Abelian link angles representing a U{1) field with residual gauge 
symmetry can then be used for the identification of monopoles. This is done, like in pure 
compact U{1) gauge field theory, by computing the U{1) gauge invariant (magnetic) fluxes 
through the surfaces of elementary 3-dimensional cubes or - what is equivalent - by searching 
for the ends of Dirac strings. Monopoles identified in this manner are generally referred to as 
DeGrand-Toussaint (DGT) monopoles We are going to localize DGT monopoles in our 
smoothed configurations put into all three gauges and to compare their possible significance. 

As stressed before, for the LG there exists an independent, gauge invariant way to localize 
monopoles in the sense of 't Hooft that can be used to test the reliability of the DGT method. 
The Higgs field introduced in the LG provides an alternative, physically more satisfactory 
possibility of monopole identification. In the continuum, lines of = \^x\ = 0, where px is 
defined as 



— Px^x ) P^ 



^E(^.")^ (8) 

\ a=l 



directly define lines of gauge fixing singularities (monopoles). For the 't Hooft- Polyakov 
monopole in the adjoint Higgs model, regions with p = of the physical Higgs field are 
identified with the centers of such monopoles. Note that this way of monopole identifi- 
cation in the pure Yang-Mills theory, too, does not require to actually perform the gauge 
transformation and an Abelian projection! It is gauge invariant by definition. 

Let us demonstrate the localization of the singularities of the LG transformation and 
how this is related to the DGT monopoles. We have plotted in Fig. ^ all lattice points 
with the local modulus squared X^(x) = J^ai^xY being less than a threshold value of 
10^^, for a generic configuration of the confined phase (t fixed, top) and for one of the 
deconfined phase {z fixed, bottom). The eigenvectors are always globally (over 12^ x 4 
lattice points) normalized to one, J2x J^ai^xY — 1- Although there are no true singularities 
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(exact zero modulus of X) located on lattice sites it can be seen that small values X^(x) 
(less than one order of magnitude below the average) mark connected regions of space on 



the lattice where one may expect to find the monopole world lines (see also Ref. |^). The 
dark lines correspond to DGT monopoles that have been extracted after Abelian projection 
has been applied to the LG. These monopoles are shown for a confinement configuration 
in Fig. ^ by arrows, where non-vertical and non-horizontal directions encode world lines 
leaving (entering) each {x, ?/)-miniplot along the z or t direction. They are located almost 
completely inside the shaded regions having X^(x) < 10~^. 

Due to the smoothness property of the gauge transformations dictated by the local 
field X^aa it has been conjectured in Ref. and demonstrated for normal Monte Carlo 
configurations that monopole currents are more abundant in the LG compared with MAG, 
as far as DGT monopoles are concerned. For RG smoothed gauge field configurations, we 
confirm this conjecture to be true as well. 

A small local modulus squared X^(x) = J^ai-^x^ expresses the obstruction to unwind 
the gauge field (to the unitary gauge) near the site x, while the local contribution to the Ritz 
functional J2a ^xi~^'^-^)x expresses the disordering of the gauge field. Here we concentrate 
on the Higgs modulus squared because this is directly related to the monopole localization. 
In Fig. H the upper plot shows the average value of provided, site by site, the local 
action SgUe falls into one of the bins on the abscissa. This plot refers to (3 = 1.4. It lets us 
expect that only points with SsUe > 4:a~'^ have a chance to meet the condition X^ < 10^^ 
mentioned before with non- vanishing probabiltity. In Fig. ^the lower plot (also for f3 = lA) 
shows the average value of X^ similarly depending on the local topological density q falling 
into one of the bins. According to this, only points with |g| > 0.01 have a chance to meet 
the X^ condition mentioned before with some probability. 

In Fig. ^ we show the probability distribution for encountering a certain X^ for the two 
cases, having no DGT monopole in a 3D box (top) and having one or more monopoles present 
(bottom), both in the confinement and in the deconfinement phase. To the (normalized) 
histograms, values of X^ found at the eight corners of each (occupied or empty) 3-D cube 
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have contributed. The histograms are averages over all configurations. It is obvious that 
in the case of a DGT monopole the distribution on the corners is strongly concentrated 
near zero. Roughly 50 % of the neighboring lattice sites fall below the threshold mentioned 
above. 

Finally, we should add one remark on the problem of potential Gribov copies. We have 
checked that the DGT monopoles, detected on RG smoothed configurations after trans- 
forming to MAG and subsequent Abelian projection, do basically not change their positions 
if random gauge transformations are applied before relaxation to MAG. From several RG 
smoothed configurations we have produced 10 different randomly gauged copies and have 
analyzed the relative shifts of individual monopole trajectories that emerge after MAG, 
Abelian projection and DGT construction. We found that more than 60 % of the dual links 
occupied by monopole trajectories will stay unchanged in position. About 30 % of the links 
will be shifted by one lattice spacing, shifts of two and more units were very rare. This gives 
reason to believe that RG smoothed fields are smooth enough such that the MAG gauge 
and Abelian projection will identify singularities without ambiguity. Notice that this was 
not the case without smoothing [§,0. Concerning the LG method, the eigenvalue estima- 
tors slightly differ numerically according to the 5 different random starting vectors that we 
considered, the regions characterized by a Higgs modulus below the threshold do practically 
not differ between the relaxed vectors. We interprete our observations mentioned before as 
an a posteriori justification for the use of the DGT method in the case of LG, too. In the 
following we will exclusively use the DGT prescription to extract LG monopoles. 



IV. PHYSICAL RESULTS ON ABELIAN MONOPOLES 



Our configurations are generated on a 12^ x 4 lattice at various /3-values with the perfect 



action tested for consistency in Ref. |3y]. In that paper we have found a critical Pc = 
1.545(10) of the deconfinement transition for A^,- = 4. According to the second order of the 
transition, we have determined it from the intersection of the Binder cumulant of Polyakov 
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lines (unblocked and blocked ones) for different spatial volumes. The present values of 
(3 = 1.4, 1.5, 1.6, 1.7 and 1.8 enclose the critical point. The smoothing is performed as 
described in Ref . . We have extracted the monopole degrees of freedom in the respective 
gauge by projection onto the diagonal part of the SU{2) gauge field. For this Abelian gauge 
field, we have then detected the DGT monopoles. 



A. Density of monopoles in various gauges 

In Fig. ^ (top) we show the temperature dependence of the average number of monopole 
currents (Mg + Mj) (total length of all monopole trajectories for a lattice of size 12^ x 4), 
for the different gauges. If monopoles are important for confinement, one would expect to 
see a decrease accompanying the phase transition. In fact, this is only the case for the LG 
and, less drastically, for the MAG. For all temperatures the monopole number in LG is 
bigger than that in MAG as expected in Ref. [^. The monopole number in PG is almost 
independent of the temperature! 

It is interesting to ask for the anisotropy ratio as a function of the temperature, 

shown in Fig. |^ (bottom). At zero temperature this ratio should be equal to one for isotropic 
gauges like MAG and LG. But it is expected to differ from one at higher temperatures. In 
our first study of RG smoothing for MAG, we observed a drop towards the deconfined 



phase much more pronounced than for unsmoothed, Monte Carlo configuration. This is due 
to the strong isotropic noise which is present in the magnetic currents before smoothing. 

We find the ratio smaller than one at our different finite temperatures, and it drops to- 
wards and across the deconfining transition. In the case of MAG, beyond the deconfinement 
temperature the anisotropy ratio is approaching zero. For LG the ratio is closer to one in 
the confinement phase and drops more abruptly crossing the transition. It decreases slowly 
for T > Tc- We have observed that also the LG monopole network decays (similar to the 
MAG monopoles [§n|)) i^^^ ^ temporally closed world lines which are, however, not 
perfectly static as in the case of MAG. In view of this, both the MAG and LG data suggest 
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the conclusion that in the deconfinement phase the corresponding Abehan monopoles be- 
come preferably static objects. What becomes static in the case of the LG are the extended 
"channels" of small Higgs modulus (see Fig. ^). For MAG, this observation has been first 
discussed in Ref. 

For the PG the anisotropy ratio of the corresponding monopoles does not depend on 
temperature at all. Moreover it is far from isotropy, with ^^jj^ ~ 0.05 (almost no spacelike 
magnetic currents) even in the confinement phase. Thus, the magnetic currents observed in 
the PG are hardly related to the confinement phenomenon. 

For clarity we show in Fig. |^ the average density of timelike magnetic currents per 3D 
hyperplane orthogonal to its direction, as a function of the temperature in units of A|. The 
upper data points represent the LG monopoles, the lower the MAG monopoles. In both of 
the gauges the physical density of the (almost static) timelike monopoles starts rising as a 
function of the temperature for T > 1.2 Tc. It has been speculated that the deconfinement 
phase is characterized by a bosonic gas of 't Hooft-Polyakov monopoles with a mass linearly 
rising with the temperature 

B. Local properties of monopoles in various gauges 

By definition, DGT monopole currents are localized along the links of the dual lattice 
which are related to 3D cubes of the original lattice. They test the properties of the non- 
Abelian gauge field most locally. We have defined in Sec. II local actions and topological 
charge densities. Now we want to discuss the different gauges from the point of view of how 
their respective monopole trajectories probe the vacuum, i.e., whether these monopoles are 
equipped with particular, gauge independent properties. We expect this picture to become 
clearer when, as in our case, UV fluctuations are removed from the gauge field configurations. 

First let us have a look how the average monopole occupation number on the dual links 
near to a site depends on the local action and charge, in analogy to Fig. ^ Fig. |^ shows 
this for the three gauges. In this respect, no drastic difference between the gauges is seen. 



15 



Again the occupation numbers of LG monopoles are bigger than in the case of MAG, those 
for PG somewhat lower. No distinction is made here with respect to the direction (spatial 
vs. temporal) of the magnetic currents. 

In Fig. |10| we show the distribution of s^-cube, separately for cubes free of monopoles 
(top) and for cubes distinguished by monopoles detected in the LG (bottom), for the lowest 
and the highest (3-value. This local action contains all loops contributing to S localized 
on the 3D cube under consideration. We see that the distributions for the two cases are 
different. Mean and variance in the case of presence of monopoles are almost a factor of two 
larger, indicating that on average significantly more action is picked up per unit of length 
along an actual monopole trajectory than would be picked up along a random walk in the 
given gauge field background. 

The situation is very similar for the MAG case, the difference between the action dis- 
tributions for monopole-free and monopole-carrying cubes is somewhat more pronounced. 
Actually, these distributions are insensitive to the direction of the monopole current only 
in the confinement phase. In the deconfined phase there are important differences. For 
LG, monopoles occupying a cube dual to a spacelike link (which are still frequent in the 
deconfined phase, in contrast to MAG!) are accompanied by a probability distribution of 
local action which is similar to that of an empty cube, and for timelike ones the action is 
not as high as in the MAG case. This explains the higher multiplicity of LG monopoles and 
the lack of suppression of spacelike "detours" for world lines closing in timelike direction. 
For cubes along the trajectories of PG monopoles the action distribution is practically the 
same as for the empty cubes. This holds in both phases. Looking back to section III.C, we 
conclude that among the gauge invariant quantities X'^ is the observable, whose distribution 
reflects the absence or presence of monopoles most clearly. This holds independently of the 
direction of the monopole current in both phases (Fig. |^). 
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C. Excess action and charge 



In a more concise way than by the distributions of local action just discussed we can 
characterize the different types of monopoles and the change with (3 by an average excess 
action of monopoles defined as 

Q ^ 'S'juonopolc 'S'no monopolc ^ /„-, 

"^ex = , \y} 

^ "J no monopolc ^ 

where S'monopoie is again the action localized on a three-dimensional cube mentioned before. 
Replacing the action in the above expression by the modulus of the topological charge density 
according to the Liischer method we can define the charge excess gex compared with the noise 
of \q{x)\ anywhere else in the vacuum. For details of the definition of the local operators 
see Ref. |20|- Fig. |lT] shows that already just below the excess action for MAG and LG 
monopoles is clearly above one (indicating an excess of action of more than a factor of two 
compared to the bulk average) and rises across the transition. Around Tc, the excess charge 
is even a few times as big for MAG and LG monopoles. These results are more pronounced 
compared with a T = study (using Wilson action, without cooling or smoothing) and 
emphasize the dynamical importance of these two types of monopoles at the deconfining 
transition. 



V. DOES ABELIAN DOMINANCE OF THE STRING TENSION HOLD IN ALL 

GAUGES? 

The Abelian dominance of the string tension in the case of MAG was a strong argu- 
ment for choosing this gauge to look deeper for the role of monopoles. Finally, this has given 
rise to the construction of infrared effective actions for QCD [|14|. Originally, of course, the 
concept of Abelian dominance was referring to generic Monte Carlo configurations. As a 
by-product of our study on monopoles in various gauges, we present here results concerning 
this question for smoothed configurations. 
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RG smoothing preserves the non-Abehan string tension within 7 % and removes the 



Coulombic part of the heavy quark force at small distance. In Fig. ^ we display the 
static quark-antiquark potential as obtained from Polyakov line correlators, for the smoothed 
SU (2) fields and, in the case of LG and MAG, after Abelian projection. The corresponding 
values for the string tensions are found in Table The Abelian string tension of the MAG 
is about 5 % less than that calculated for the non-Abelian SU{2) field. The Abelian string 
tension of LG is again somewhat smaller than for the MAG. In qualitative respect, we still 
find Abelian dominance for smoothed lattice fields, independent of the gauge chosen for 
projection. Notice, however, that in the PG (where the temporal links f/(x,2:4),4 become 
diagonal and X4 independent) the non-Abelian string tension extracted from Polyakov line 
correlators is trivially identical with the Abelian one. 



VI. CONCLUSION 

We have collected additional evidence that the RG smoothing technique (with an ap- 
proximate classically perfect action) is a powerful tool, here used to investigate semiclassical 
aspects of the monopole structure of the Yang-Mills vacuum. We have presented further 
results concerning the recently proposed, physically better motivated and gauge invariant 
method to identify monopoles on the lattice, the so-called Laplacian gauge (LG). We con- 
trasted the corresponding monopoles with other types of monopole trajectories which usually 
are obtained in a gauge dependent way. What parallels exist and how different the monopole 
content can be if different gauges are used to localize them has been examplified by consid- 
ering three gauges (MAG, LG and PG). In a first step we showed that the gauge invariant 
way to localize monopoles is consistent with the DeGrand-Toussaint (DGT) method being 
applied to Abelian projected configurations after LG fixing. This correspondence becomes 
clearer if smoothed Monte Carlo configurations are analyzed. This justifies to take over 
the DGT method to study the monopole content of configurations also in the LG, at least 
after RG smoothing. Accepting this technique, we also found strong correlations between 
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positions of MAG and LG monopoles. 

While MAG and LG monopoles (as extracted by the DGT method) behave similar at 
the deconfining phase transition concerning density and anisotropy, monopoles identified 
in the Polyakov gauge apparently lack a dynamical relevance for confinement and at the 
deconfinement transition. In a second step we then analyzed trajectories of monopoles 
identified in various gauges and found that monopoles (in all three gauges studied) appear 
preferably in regions which are characterized by enhanced action and topological charge 
density. The reverse is not universally true: only monopoles detected in LG and MAG 
are characterized by significantly different local distributions of action. We have shown 
earlier that strong gauge fields are locally (anti) self dual. This is further, circumstantial 
evidence that the monopoles are in fact dyons in the confinement phase. As a by-product 
of our study we found that almost the complete string tension can be recovered from the 
Abelian projected field corresponding to various Abelian gauges, including LG and PG. For 
the latter this is trivially true, as long as we measure the string tension by means of the 
Polyakov line correlator. 

We have shown that for smoothed gauge field background monopole trajectories carry 
an excess action of about twice the bulk average of local action. Similarly, monopoles also 
carry an even bigger excess topological charge. In the confinement phase this observation 
applies to monopoles independent of which gauge has been chosen for identifying them, but 
PG monopoles are very different in the deconfinement phase. We conclude that the DGT 
monopoles related to MAG and LG behave similar physically and can be semiclassically 
interpreted as physical objects which carry considerable action and topological charge in the 
vicinity of the deconfinement transition. 

This work was supported in part by FWF under project P11456. One of us (M. M.- 
P.) acknowledges support by the European TMR network Phase Transitions in Hot Matter 
under contract FMRX-CT97-0122. E.-M. I. wishes to thank T. Suzuki and M. I. Polikarpov 
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TABLES 



w{i.j) 


./ = 1 


./ = 2 


./ = 3 


J = 1 


i = 1 (plaquettes) 


1.115504 


-.5424815 


.1845878 


-.01197482 


i = 2 (6-link loops) 


-.01443798 


.1386238 


-.07551325 


.01579434 



TABLE I. Weight coefficients of the simphfied fixed-point action. 
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variants of pure SU{2) theory a 

non-Abelian, without smoothing 0.22(2) 

non-Abehan, with RG smoothing 0.205(2) 

RG smoothed, projected to Abehan from MAG 0.189(1) 

RG smoothed, projected to Abelian from LG 0.180(2) 

RG smoothed, projected to Abehan from PG 0.205(2) 
TABLE II. String tensions extracted from Polyakov hne correlators at /? = 1.4. 
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FIG. 1. Abelianicity as a function of f3 for MAG, LG and PG. 
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FIG. 2. Average over 100 independent configurations of the lowest eigenvalue of the adjoint 
lattice Laplacian on the 12^ x 4 lattice as a function of (3 or temperature. 
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FIG. 3. Regions with local Hiegs norm 

Y^ai^x)^ < 10"^ (marked by dots) of the lowest 
eigenmode of the adjoint lattice Laplacian (normalized over the 12'^ x 4 lattice). Monopole 
loops (shown as lines, for clarity shown also separately in the right plots) have been identified 
following DGT after Abelian projection from the Laplacian gauge. The upper configuration is 
taken from the confinement ensemble for a fixed time slice, the lower from the deconfined phase 
for z fixed. 
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Beta = 1.4, Q = 1.79 
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FIG. 4. Contour plots in the (x, y)-plane of the local Higgs norm (shown only for 
Y.ok^x)'^ < 10"^) in the z = 1,..,6 half of the lattice for all time slices t = 1,..,4. Magnetic 
currents of DGT monopoles as obtained by Abelian projection from LG are symbolized by arrows 
with non- vertical or non-horizontal directions indicating the z and t directions. 
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FIG. 5. Average of the local Higgs norm X'^{x) = J2ai-^x)'^ lattice points having a value of 
local action SgUe within the same respective bin (top) . The same for lattice points binned according 
to topological charge density q (bottom). 
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FIG. 6. Probability distribution of the local Higgs norm in the absence (top) and in 
presence (bottom) of a DGT monopole obtained in the LG. f3 = 1.4 (left) refers to the confinement, 
/3 = 1.8 (right) to the deconfinement phase. 
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FIG. 7. Average total monopole length on the configurations (top) and the space-time asym- 
metry of magnetic current densities (bottom) shown for three different gauges as a function of (3. 
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FIG. 8. Average spatial densities of (timelike) monopole and antimonopole currents shown for 
MAG and LG as a function of physical temperature. 
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having a value of local action SsUe within the same respective bin (top). The same for lattice 
points binned according to topological charge density q (bottom). Both plots refer to f3 = 1.4 
(confinement phase). 
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FIG. 10. Probability distribution of cube action s^^cube in the absence (top) and in presence 
(bottom) of a DGT monopole obtained in the LG. P = 1.4 (left) refers to the confinement, /? = 1.8 
(right) to the deconfinement phase. 
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FIG. 11. Excess action (top) and topological charge (bottom) shown as function of physical 
temperature. 
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FIG. 12. Heavy quark potential evaluated after smoothing from non-Abelian SU{2) Polyakov 
line correlators at /? = 1.4 (confinement) and from the Abehan Polyakov line correlators in Abelian 
projections from MAG and LG, respectively. For PG there is no difference with the non-Abelian 
potential. 



36 



